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SUMMARY 


The present report presents the theoretical formulation, analysis 
and results necessary to analyze quadrupole noise generated from a loaded, 
isolated, subsonic rotor due to its interaction with an inflow distortion 
or inlet turbulence. Programs to calculate these interactions are given 
in an appendix. These programs also yield dipole noise contributions 
from these two mechanisms. The ratio of quadrupole to dipole noise is 
largely a function of the axial flow Mach number, wheel tip Mach number, 
rotor solidity, total pressure ratio across the rotor. It is relatively 
independent of the specific form of the inflow distortion or inlet 
turbulence. Comparisons with experimental data only succeed in predicting 
gross levels at a given speed and fail to predict the variation of noise 
at fixed speed with flow and pressure ratio. Looking at the separate 
variations of the theoretical predictions of dipole and quadrupole noise, 
however, one notices that if these levels were closer together, a suitable 
composite of them would yield the trends of the experimental data. This 
leads to a suggestion that the calculation may be overestimating the ratio 
of dipole to quadrupole noise by overestimating the dipole noise. Likely 
sources of this overestimate are suggested along with corrective 
procedures . 



INTRODUCTION 


The purpose of the present study is to quantify a quadrupole noise 
source which was first proposed in (2) as a likely contributor to fan 
noise. The source arises from the fluctuating Reynolds stresses intro- 
duced by a combination of two unsteady velocities, one being due to the 
potential flow field of the rotor and the other due to a solenoidal velo- 
city field such as an inlet distortion or inlet turbulence. 

The study proceeds as follows. Since interest is in fan/compressor 
noise where substantial duct axial velocities are involved, a pertinent 
Lighthill equation including uniform axial flow is derived first for the 
study of quadrupole noise. Secondly the potential flow fields of a sub- 
sonic rotor are derived as a function of the pressure ratio across it, 
the wheel tip' and axial flow Mach numbers. Thirdly a systematic approach 
to the acoustic problem of estimating the noise from this source is worked 
out based on repeated use of Fourier exponential transforms. 

We next present calculations under several constraints of the 
relative magnitudes of the quadrupole versus the dipole noise source. 
Comparisons are also carried out with some available data on the in- 
fluence of pressure ratio on fan noise. In an appendix, computer pro- 
grams for calculations of quadrupole and dipole noise due to both inlet 
distortion and inlet turbulence are given. 

Ivan H. Edelfelt of the General Electric Research and Development 
Center helped considerably in programming these calculations. Thomas F. 
Gelder and Marvin E. Goldstein of the NASA-Lewis Research Center provided 
several helpful discussions. 



DEVELOPMENT OF PERTINENT LIGHTHILL EQUATION FOR BLADE LOADING NOISE 


The full Lighthill equation for aerodynamic noise for an inviscid 
gas, in the absence of mass, energy or force sources (and assuming the 
isentropic relation p = a^ p) may be written as: 
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Now let u^ = U 6-|^- + uj where U is uniform and steady and 1 corresponds 
to the axial direction. Then we may show rigorously that (1) reduces to: 
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Equation (2) is the desired form of the Lighthill equation for blade 
loading noise from an isolated rotor. The form of p u^ u j , correct to 
second order, would be (pQ u| uj). For an axial flow device, two im- 
portant sources of velocity fluctuations are solenoidal velocity fluctua- 
tions associated with inlet distortion or inlet turbulence and secondly 
irrotational velocity fluctuations associated with the flow field of a 
loaded rotor. Thus we may write u^ = u^g + u-jp where u|g denotes the 
solenoidal velocity fluctuation associated with inlet distortion or turbu- 
lence and Up that associated with the rotor potential flow field. The 
calculation^of the quadrupole noise field due to interaction between the 
rotor potential flow field and an inlet distortion or turbulence field 
thus boils down to solving (2) with a source term of type, 


ax. ax. ^^0 *^ip 


The factor of 2 arises from the symmetry of the quadrupole term with 
respect to interchange of i, j. 

DETERMINATION OF FLUCTUATING VELOCITIES DUE TO THE ROTOR POTENTIAL FLOW 
FIELD 


For the determination of the fluctuating velocity field of the rotor, 
the most convenient quantity is the (nondimensional ) lift coefficient of 
the steady lift exerted by the rotor blades on the fluid. However the 
quantities most conveniently available in practical terms are the pres- 
sure ratio across the rotor, its wheel tip and axial flow Mach numbers, 
etc. So we first evaluate C|_ in terms of these quantities. 

The work per Ibm of the flow is 
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is gives the desired expression for C|^ as: 
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To calculate the u', v' (x and y components of velocity) requires a 
linearized, compressible analysis of the flow field of the isolated rotor. 


Let Wy. = a Mf,. We first find the solution for velocity components 
u', v', parallel to the x' - y' coordinate system of Figure 1, due to 
equally spaced concentrated unit forces at the origin and its corre- 
sponding points as shown in Figure 1. (Note that the blade exerts a 
force on the fluid equal and opposite to the force by the fluid on the 
blade.) We use a frame of reference fixed w.r.t. the translating blade 
row so that we have a steady state problem. We have to consider the 
effect of a sum of forces: 
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where 1_ denotes a unit force vector, and 6 stands for the Dirac delta 
function. By using a result on page 68 of (1) concerning the sum of an 
infinite row of equally spaced delta functions, clearly the above is equal 
to: 
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The linearized equations of motion and continuity are: 
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Eliminating p' from (lb, c) we derive that: 
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and since u', v' and the delta function term vanish far from the blade row 

(2) 
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[Kutta Joukowski Law] 

Next we eliminate the force terms in (lb, c) and assuming an isentropic 
relation between p' and p' one obtains: 

[1 - cos^(a^)] + 1^ [1 - sinMa^)] 

= Mj sin(a^)cos(a^)[|^+ 1^] 

(modified continuity equation) 


(3) 


Using (2), (3) single equations for u', v' may be obtained which may 
be solved by requiring that u’, v' vanish as x' ^ ± We omit the details 
and give the result. 

u' for x' < 0 = 
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Similarly: 

v' for x' < 0 = 
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*sgn(n) =■ 1 if n > 0, = -1 if n < 0 and = 0 if n = 0. 
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When viewed in a frame of reference fixed with respect to the casing of 
the machine, these u', v' velocities given by (4a), (4b) appear as: 
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v' for X < 0 = 
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THE AXIAL VELOCITY OF THE INLET DISTORTION 


Far upstream from the rotor, the inlet distortion appears as regions 
of axial velocity defect superposed on a uniform inlet velocity. 

Some assumption regarding the form of the inlet distortion is 
necessary in order to analyze its effect. The form assumed is sketched 
in Figure 2. 

Let the axial velocity, in Mach number, associated with the uniform 
total pressure region be M^. The axial velocity defect AM, associated 
with the inlet distortion may be deduced as follows from the condition 
that far upstream there is no variation of static pressure. 
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(y is the specific heat ratio of the gas) 
Now 
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Thus one may solve for AM, in (6) as: 
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The form of Eq. (8) indicates that Ap./p. and M are subject to a 
restriction that ^ 


(1 



(1 + %) > 1 


(9) 


(Other than the above restriction, M and Ap./p. may be specified 
arbitrarily.) 


The restriction of (9) follows from the assumption that the dis- 
tortion consists of regions of constant static pressure and varying axial 
velocity. For such a distortion, (9) simply expresses the requirement 
that the static pressure not exceed the total pressure anywhere. For 
small Ap./p. , and small Mf, (9) may be expressed approximately as 

L L a 





(9) 


The interaction of the inlet distortion with the isolated rotor may 
be analyzed (as sketched in Figure- 2) as the interaction of a shear wave, 
convecting along the fan inlet duct at an axial Mach number M , with an - 
isolated rotor. . ■ • 

Let the inlet distortion be Fourier analyzed in an x-y coordinate 
system at the mean radius of analysis "a" as follows. The basic 
periodicity of the distortion in the y-di recti on is "Zira'VS (for an S 
lobed distortion). Let the distorted velocity profile be expressed as: 
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THEORY OF QUADRUPOLE NOISE GENERATION BY INLET DISTORTION 
C a W M M 
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2B (Sometimes B is also used to denote 
the number of rotor blades and the 
context of its use makes this fairly 
clear.) 
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Only an axial quadrupole of type p u u and a transverse one of type p 
V u where u , v are the axial aHd tangential velocity distortions 
aiso^iated wi?h tRe rotor potential flow field and u the axial velocity 
distortion associated with the inlet distortion contributes to the noise 
(at blade passing frequency and its multiples). 

Exploiting the three facts that (a) there is no interest in the 
zeroth harmonic of blade passing frequency (which is a dc field) (b) that 
the rotor is subsonic (M^ < 1) and (c) we are only interested in pro- 
pagating noise fields, we may show that the noise generation for a 
particular harmonic "n" (^ 1) of blade passing frequency noise is governed 
by: 


* "a P - '“P 


KiA' * j B') I q 

1 ■ 1 


exp(-a^|x| )exp(j 6^ x)exp(j B^(l - 


- 3^ 

+ 3x 3y ^0 


C 


i 


I C: exp(-a^|x|)exp(j 6^ x)exp(j 3^(1 - ■^)y) 
= 1 


7 


/ 



(where C ! = a C . ) 

Clearly dep^nd^ on y, t (given i) as 

exp[j B^(l - ^)y]exp[j ' 


Thus we deduce that the axial part of p depends on x as: 
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To solve (2) introduce the Fourier exponential transform P(z) of p by: 
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By the method of residues, the downstream/upstream waves have amplitudes 
given by: 



( 5 ) 
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Independent waves (above cut-off) for each n will be produced over a 
range of i given by: 
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For each of these i, amplitudes of upstream and downstream waves are 
computed from (5). It is demonstrable that the rms values of the axial 
components of the intensity are (power/cross sectional area of annulus): 
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A = [M2 - (1 - m2)(1 - 

This completes the theory of inlet distortion quadrupole noise generation. 


NOISE GENERATION BY INLET TURBULENCE 

We represent the turbulent velocity components by: 


(Uj, Vj, Wjl 


u,v,w ' ' 


all k 

adopting the usual Fourier Stieltjes form where 


dZ^(k' )dZ^(k‘ ' ) = 6(k' - k" )o^y(k’ )d^k' d^k’ 


and 'l’uy(k) is the usual spectral density tensor. 
The noise generation problem is given by: 
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A typical term of the RHS of (7) may be written as: 
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( 8 ) 


-pQ Ip- U±A + sgn(n)j B)dZu(...)} 

■Oo ifV j sgn(n)B)dZ^(...) 

• C dZ„(...)l ; p^ l^C dZ,(...) 
where (...) stands for: 

exp[-a^lxj]exp[j((S^ + k^)x] 

exp[j(B + k )y]exp[j[6„ M - M k la t] (9) 

n y n L d X u 

A word of explanation is in order here concerning the model of 
turbulence employed in the analysis. The model is the same as was 
employed in (3). A three dimensional model of turbulence is employed 
but the spectrum functions are integrated over the third coordinate 
(z) so that nowhere does the z dependence appear explicitly. 

Substituting for the y dependence from (9), the source term (8) is 
effectively: 

■^0 Ix^ 

-Po I7 j(8n l<y){(±A + j sgn(n)B)dZ^ ± C dZ^)(...) 

± (6n + Pq C dZ^(...) (10) 

As before the Fourier exponential transform of the term in (10) may be 
written as: 
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where the definitions of G (z) and G (z) are obvious. 
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As before the downstream and upstream waves will have form: 

(i) if (B„ Mx. - k^ M, ) > 0, then with 
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and the rest is as in (13). 

Let A = [(B* M^. - k M )2 - (1 - m2)(b + k 
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f = 13, k^Mj (15b) 

Then the acoustic energy produced downstream and upstream by such waves 
is given by taking the mean square of (13) and multiplying the resulting 
quantity by: 
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The remaining procedure is more or less mechanical. Let us say that'the 
interest is in the acoustic energy between wave numbers lying between x 
and X + Ax- (Wave number = frequency in radians/sec. speed of sound). 
This determines that for each n from -«> to <=°, a range of k lying, between 
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is of interest. This determines (for each n) a specific value of k . A 
range of from ^ 
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is to be considered (depending on whether x is positive or negative). The 
quantity 


GMz_^)< 1> + 2G (z_,)G (z^)$ + GMz_^)iI> 
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times the factor in (16) is to be integrated over the range of k indi- 
cated in (18), given x> 6 (and hence k by (17)) for each nin ^ange 
to «>. For given physical^ave number o^ interest both energy corresponding 
to +x and -x need to be considered. By adding up the infinite series of 
contributions from n = to n = «> for +x and -y, taking note of the fact 
that the width in k is (Ay/M,), etc., the required spectral density 
dl-/dy may be deduced. $ , $ and <I> should first be integrated over 
k from -<» to <», in keepiHy wi'iK two dYI¥iensional or plane nature of present 
analysis. Actually it turns out to be more pertinent to think in terms of 
X dl±/dy as this quantity has the units of power and is related directly 
to measurements obtained by constant percent bandwidth or constant octave 
filters. If we assume a specific form of the longitudinal velocity corre- 
lation function for the turbulence of type exp(-r/L), $ » ® ^vv 

can be easily written down for given turbulence intensii^y. ^ 


DISCUSSION OF RESULTS 

Previously developed analyses (3, 4) of dipole noise were employed 
to estimate the dipole contribution to inlet distortion and inlet turbu- 
lence noise. 

We show in Figures 3(a), (b), (c) first the nohdimensional constant 
percent filter spectra due to inlet turbulence rotor interaction noise. 

We study a fixed operating line characterized by a steady lift coefficient 
of unity associated with a rotor whose solidity is unity. In turbu- 
lence noise calculations', one parameter that enters the calculation is 
the ratio of the integral length scale of turbulence to the blade spacing 
designated herein as (L/D). We have independently verified that the 
specific value of the ratio (L/D) does not materially affect the ratio of 
dipole/quadrupole noise contribution. For rotors with axial Mach numbers 
of order two-thirds the wheel tip Mach number, at wheel tip Mach numbers 
of order 0.8, the quadrupole noise contribution does indeed begin to 
exceed the dipole contribution especially at the higher frequencies. The 
low frequency end of the spectrum (frequencies less than half. the blade 
passing frequency) is still dominated by dipole noise. Similarly calcu- 
lations were carried out Tor a rotor with (Ma/M-t) = 0.5 and C|_ = 0.5 and 
are shown in Figures 4(a), (b) and (c). Similar trends are evident insofar 
as the fact that the ratio of quadrupole/dipole contribution increases 



with M|. and also, at a given increases with frequency. However one 
finds that the ratio of quadrupole/dipole noise, at given M^-, increases 
with the ratio To bring this out in somewhat sharper focus, a 

calculation is shown in Figure 4(d) of a rotor with (Mg/M^) being unity, 
at the same relative Mach number and lift coefficient, etc., as the rotor 
in Figure 4(c). The ratio of quadrupole/dipole noise is much greater in 
Figure 4(d) than in Figure 4(c). In Figures 5(a), (b), (c) and (d) are 
shown calculations of noise spectra at a fixed axial Madi number and total 
pressure ratio but with varying tip speed. Again as might be expected, 
higher tip speeds (and associated higher frequencies) do increase the 
relative ratio of quadrupole/dipole noise. However no evidence is 
available in these results to support the notion that fan noise could be 
minimized for constant work by choice of high tip speed and low loading. 
For other specific designs over some limited tip speed ranges, some such 
result may be true but it is certainly not a generally valid principle. 
Indeed the results of Figure 5 support the notion that for subsonic fans, 
to minimize noise for a given flow and total pressure ratio one should 
employ high blade loading and low tip speed. This conclusion applies of 
course only to the physical generation process and not to aspects such as 
effectiveness of treatment or perceived noise considerations. 

In Figures 6 and 7, calculations corresponding to Figure 5 and Figure 
3 are shown for 1st and 2nd harmonic pure tone (blade passing frequency) 
noise generated by an inlet distortion - rotor interaction. In both cases 
a maximum total pressure defect of 2% of the inlet total pressure was 
assumed. A four lobed distortion is assumed to impinge on a 38 bladed 
rotor. Other details are given in Figures 6, 7. A fairly smooth profile 
of the distortion as shown in Figures 6, 7 was assumed. The results of 
Figures 6, 7 are similar to those of Figures 5, 3 in that the ratio of 
quadrupole/dipole noise is greater for second harmonic rather than first 
harmonic noise and also increases with tip Mach number when a fixed 
operating line study is considered (ratio of (Mg/M^-) fixed, fixed C|_). 

Finally an attempt was made to compare the results of the current 
analysis with experimental data from (5). The data given in (5) pertains 
to 50 Hz bandwidth forward radiated power at the blade passing frequency 
measured at several speeds and several pressure ratios at each speed. 

The data for two rotors (designated rotor 1 and rotor 2 in (5)) were 
employed. 

Figure 8 indicates the results for rotor 1. This case has also been 
analyzed by the authors of (6) and, as assumed by them, a four lobed dis- 
tortion with maximum velocity defect equal to \% of the axial velocity is 
assumed to impinge on the rotor. The shape of the distortion was assumed 
to be a triangular pulse of width 10% of the extent of the lobe. This 
width corresponds roughly to the width of each of four struts placed up- 
stream of the rotor in the experiments reported in (6). 

Shown in Figure 8 are the actual measured data at 50, 60, 70 and 80% 
speeds and predictions by the current analysis using the results of the 
current study and those of (4). The dipole and quadrupole contributions 
are shown separately (as predicted by the theory). The dipole noise 
levels predicted by the theory are roughly in the ballpark of the results 
observed in J5) by Gelder and Soltis but in one major respect the theory 
fails completely. The present calculations employ the two dimensional 
Sears gust formula for estimation of unsteady forces needed to obtain 
the dipole noise. According to these results, the dipole noise greatly 
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exceeds the quadrupole noise (by at least 10 dB) over the operating range 
of rotor 1 of (5). Since dipole noise is independent of blade loading 
and is proportional to the rotor relative velocity (which decreasea as 
the loading is raised at constant wheel tip speed), the inlet noise is 
always predicted to decrease with increasing pressure ratio (and 
associated lower weight flow and rotor relative velocities). .The data of 
(5) of course fail to show this trend at all. 

In case of rotor 2, an attempt was made to predict the noise (in a 
50 Hz bandwidth) with the aid of the turbulence noise prediction pro- 
cedures. Measurements of scale and intensity of turbulence are pretty 
scarce but a good set of measurements for one flow condition and one inlet 
is available in (9) and these were crudely scaled to the conditions of the 
inlet in (5). A 4% intensity of turbulence effective over the outer 10% 
of the rotor blade span and an (L/D) of about 0.5 was assumed for these 
calculations. Again one notices a fair ability to account for the gross 
noise levels by the inlet turbulence dipole noise theoretical estimates 
but a total failure to predict the variation with pressure ratio and flow 
at fixed speed largely due to the extreme predominance of dipole noise 
according to the theoretical estimates. 


To put these results in perspective, we start by remarking that the 
quadrupole noise mechanism pointed out in (2) is a second order effect 
while dipole noise is a first order effect. Thus, with a dipole mechanism, 
the acoustic pressures scale roughly as (po Wp u') where Pq = mean fluid 
density, Wp = steady relative velocity through the rotor and u' the 
fluctuating solenoidal velocity (whether due to distortion or turbulence). 
With a quadrupole noise mechanism, they scale as (po Up u') where Up is 
the potential flow field induced by the rotor and generally Up << Wp. 

Indeed the very procedure adopted to estimate Up is a small perturbation 
calculation assuming that Up and Vp « Wp. 


Again the full quadrupole source term » ^ — (p„ u. u.) may be 

expanded as (0^ + 2u.j + ( 9 ^)( 9 ^)) where 0 is the divergence of the 

3 

velocity. So far as the inflow distortion or inflow turbulence is con- 
cerned, these velocity fields have zero divergence. Even the potential 
flow field of the rotor does not have a significant divergence at modest 
axial and wheel tip Mach numbers. For Mg = M+, for example, the divergence 
of the rotor velocity field is M|(0 + 2(9u/9y)) and thus of second order 
in the axial Mach number. Hence only the term 


9u. 9u. 


9xh 


(p u. u.) may be expected to contribute significantly. 

^ * V 


Finally the quadrupole source is an extended source subject to extensive 
phase cancellation while the dipole source is generally assumed to be 
compact (as is the case in the present model). In this regard it is 
interesting to note that the authors of (6) considered only the potential 
flow field upstream of the rotor plane as contributing to the upstream 
noise and it is possible that such an approach will involve less phase 
cancellation. 


Rao and Chu (8) concur with the results of the current calculations 
in obtaining estimates of the ratio of dipole/quadrupole noise of the 
order of 20 dB. Morfey (7) has calculated noise from the axial quadru- 
pole (pq Up Us)but it is difficult to see why he did not include the 
contribution of the x-y quadrupole (p v u ). Thus the term that Morfey 
calculates is P 



^2 ^2 

^- 7 -(u„ u ). Actually if one includes the term v„ u^) the 

0 dX p $ ’ dX dy ops 

quadrupole noise source term can be written as. 

’“s * 5, 

^0 IT ^ ‘’o S ar 


where 0 = divergence of rotor potential flow field = 

3u Sv ' 

(■r-^ + • 5 —^). For the low Mach number applications considered by Morfey 
dX dy 


(he has 0 < M. < M and M < 0.7) and especially for the Mach numbers 

^ t ““ d d 

(M ~ 0.25) where he claims the quadrupole mechanism overtakes the dipole 
mechanism, 0 p was probably negligible and thus the term he should have 
dealt with was really 

3v - g2 

■ 5 : 7 ^ • The term that he really employed viz. p„ ^^^(u u ) is 

0 dX dX 0 dX P S 

actually more or less cancelled by a part of (Pq Vp u^). Thus a 


serious question arises with regard to the correctness of his relative 
estimates of the quadrupole/dipole mechanism. 


We are still left however with the question of the total failure of 
the present results in Figures 8 , 9 to agree with the data of (5). Two 
remarks are in order here. First of all, absolute prediction is really 
much less of a goal than the prediction of relative trends. This is 
because absolute prediction requires a whole host of good inputs other 
than the acoustic data such as distortion, turbulence data which have 
only been estimated in this case. Paradoxically of course the present 
predictions have done somewhat better on absolute prediction than on the 
relative trends. We notice consistently from Figures 8 , 9 that if the 
dipole-quadrupole levels were more closely matched, the relative trends 
in the data would be much better explained. In other words, the actual 
data in Figures 8 , 9 often exhibit trends intermediate between those of 
the dipole and quadrupole noise trends considered separately. Needless 
to say, since the dipole noise is predicted to be so much higher than the 
quadrupole source, any attempt to simply add them algebraically would of 
course fail to explain the data trends. This brings us to the second 
remark. It appears that the present calculations may be considerably 
overestimating the ratio of dipole/quadrupole energy. The quadrupole 
noise calculation is, in many respects, the more exact calculation. 

The dipole noise estimate ertjploys the Sears gust formula and there are 
indeed several reasons for believing that it leads to overestimation of 
the unsteady force. Note that the Sears gust formula applies only to 
plane, incompressible flows involving isolated, flat plate airfoils at 
zero angle of attack. Firstly accounting for longitudinal or chordwise 
gusts generally leads to reduced values of gust loads (7). For inlet 
distortion noise, the reduction may be estimated as (tan(a ) - C./tt)/ 
tan(a^). Secondly suggestions have been made that the Sears function 
could be revised to accommodate real fluid effects by multiplying it by 
the steady value of (dCL/da) at finite loading to the steady value of 
(dCL/da) at zero loading. This correction would have to be determined 
experimentally but (dCi/da) generally decreases with increasing loading 
(or C|_) becoming zero Before the airfoil stalls. Similarly inclusion of 
finite aspect ratio effects or of gust obliquity would also lead to 



lowered estimates of unsteady loads (see (6) for suggestions on including 
the aspect ratio effect - the suggested procedure is very similar to the 
one concerning real fluid effects). Finally compressibility and cascade 
effects also generally tend to lower estimates of the gust loads. The 
work of Timman (10) on compressibility effects at high Mach numbers for 
isolated airfoils does predict a reduced gust load (at fixed reduced 
frequency). The source compactness assumed in the present calculations 
for dipole noise may again be a source of overestimation in view of the 
tendency of phase cancellation to occur with a distributed source. 

One final remark with regard to Morfey's study (7) is that by choice 
of a distortion 'with period equal to the blade spacing (one that is not 
likely to occur in practice), he starts, a priori, with a situation of 
low gust loads due to the very high reduced frequencies involved. Also he 
considered a rotor of unusually low tip solidity = 0.5, a model which 
again inherently favors the quadrupole mechanism. 

CONCLUSIONS 

1. A systematic theory of quadrupole noise generation in fans/compressors 
has been worked out both for inlet distortion and inlet turbulence - 
rotor noise. The calculations along with previously developed analyses 

of dipole noise due to these agencies have been coded into working pro- 
grams included in the present report. 

2. Parametric calculations with hypothetical constraints indicate that; 

a) At any fixed operating point, the ratio quadrupol e/di pole noise 
increases with frequency whether one considers inlet turbulence 
or inlet distortion to be the source of noise. 

b) Along a fixed operating Vine, i.e. for fixed (M^/M^) and for 
given C|_, the ratio of quadrupole to dipole noise increases with 
tip Mach number. 

c) At fixed Mg and given total pressure ratio, the variation of 
ratio of quadrupole to dipole noise with tip Mach number is not 
very simple. This is because with increasing tip speed and 
fixed pressure ratio, the loading (or Cj^) decreases with tip 
speed but the frequencies increase. These features are illustrated 
in Figure 5 where the ratio of quadrupole to dipole noise is 
minimum in Figure 5(c) at = .6. 

d) The ratio of quadrupole/dipol e noise is largely a function of 
Ma> M^, rotor solidity and pressure ratio (or C^) and largely 
independent of the specific form of the distortion or turbulence 
(within broad limits). 

e) The analysis shows the quadrupole noise to be independent of 
rotor solidity (at fixed Mg, M^ and total pressure ratio) but 
dipole noise increases with rotor solidity and hence the ratio 
of quadrupole to dipole noise decreases with increasing rotor 
sol idity. 

3. No general support is obtained for the idea that for given total 
pressure ratio, fan noise can be lowered by operating at high tip speed 
and low loading based on present calculations (see Figure 5). 



4. Comparisons with data from two rotors as reported in (5) reveal that 
the analysis does fairly well in predicting the gross PWL levels at given 
speeds but fails substantially in predicting the trends of the variation 
of PWL at fixed speed with flow and pressure ratio. When this deficiency 
is examined in detail, it is found that if one looks at the separate 
theoretical predictions of the trends of variation of PWL due to dipole 
and quadrupole mechanisms, the data of (5) often exhibit trends intermediate 
to those of the dipole and quadrupole predictions. However, the theory 
consistently places the dipole noise levels at values much higher than the 
quadrupole levels for the design parameters of the two rotors studied in 
(5) and hence a composite of the two theoretical predictions does not 
exhibit the trend of the real data of (5) at all. A possible inference 
from this is that the ratio of dipole/quadrupole noise is overestimated 
by the theory and it is argued that this may be chiefly due to overestimate of 
the gust loads by relying on the Sears gust formula to predict them. 
Modifications to the gust formula to allow for real fluid, chordwise gust, 
aspect ratio, spanwise variation, compressibility and cascade effects are 
all suggested refinements likely to bring the predicted levels of dipole 
and quadrupole noise much closer together. If this happens, it is possible 
that a composite of the two will result in predicted trends in much better 
agreement with the data of (5). 



APPENDIX 1 


COMPUTER PROGRAMS FOR CALCULATION OF DIPOLE AND QUADRUPOLE 
NOISE DUE TO INLET DISTORTION OR TURBULENCE 


INLET DISTORTION NOISE 


Inputs needed are: 


TPD 

IS 

EMA 

EMT 

NDP 

6AM 

RH0 

C 

A 

N 

IB 

SIGR 

TPR 

TPDS(I) 


total pressure defect, maximum value of (Ap^)/p^. 
number of lobes in distortion. 

^a 

''t 

number of points used to specify the distortion 
shape (see Figure Al), should be less than 51. 

Y 

gas density in Ibm/cft. 
speed of sound in fps. 

area of annulus being studied in square feet, 
harmonic of blade passing noise of interest, 
number of rotor blades. 


a 


^^ 02 ^^ 01 ^ 


NDP values specifying relative shape of total 
pressure distortion (max. value would be unity: 
0 < TPDS(I) £ 1). See Figure Al. 


Outputs given are:' 

1) All inputs except NDP, TPDS(I). 

2) Sound power in dB re: 10" watts in each lobe number M 
where M = N • IB - I • ( IS) . M, I are also printed. The 
dipole and quadrupole noise upstream and downstream are 
given as PUDDB, PDDDB, PUQDB, PDQDB. The sum over all 
propagating M is given as 

SPUDDB (dipole noise upstream), 

SPDDDB (dipole noise downstream) , 

SPUQDB (quadrupole noise upstream) and 

SPDQDB (quadrupole noise downstream). 



APPENDIX 1 


INLET DISTORTION NOISE 


Inputs needed are: 

TPD total pressure defect, maximum value of (Ap^)/p 

IS number of lobes in distortion. 

EMA M 

d 

EMT 

NDP number of points used to specify the distortion 

shape (see Figure Al), should be less than 51. 

GAM Y 

RH0 gas density in Ibm/cft. 

C speed of sound in fps. 

A area of annulus being studied in square feet. 

N harmonic of blade passing noise of interest. 

18 number of rotor blades. 


SIGR a 

TPR 

TPDS(I) NDP values specifying relative shape of total 

pressure distortion (max. value would be unity: 
0 5 TPDS(I) £ 1). See Figure Al. 

Outputs given are: 


1) All inputs except NDP, TPDS(I). 

-13 

2) Sound power in dB re: 10 watts in each lobe number M 

where M = N • IB - I*(IS). M, I are also printed. The 
dipole and quadrupole noise upstream and downstream are 
given as PUDDB, PDDDB, PUQDB, PDQDB. The sum over all 
propagating M is given as 

SPUDDB (dipole noise upstream), 

SPDDDB (dipole noise downstream), 

SPUQDB (quadrupole noise upstream) and 

SPDQDB (quadrupole noise downstream). 



INLET TURBULENCE NOISE 


Inputs needed are: 


EMA 

EMT 

TPR 

NELMAX 

M 

TI 

GAM 

RH0 

CS 

AI 

SIGR 


M 




needed to be investigated, 
dl. 


number of (L/D) 
number of frequencies at which f is desired, 
turbulence intensity <u'^> v U^ 

Y - 

p in Ibm/cft. 

speed of sound in fps. 

cross sectional area of annulus of interest, sq. ft 
a 


F(I)(I = 1, M): M values of frequency expressed as (f/fu) 

where fj^ = blade passing frequency. 

ELIN(I)(I = 1, NELMAX): NELMAX values of (L/D) of interest. 


Outputs are: 

1) All inputs except NELMAX. 

5 

2) For each L/D (written as EL) and F, 10 log^^^If ^ "2 P 

for both dipole and quadrupole sources as SDVDBR, SDDDBR, 
SQVDBR, SQDDBR, Also sum of dipole and quadrupole 

j T IQ 

contributions to [f ^] , upstream and downstream re: 10''^ 

watts as PVDB and PDDB. 



INLET TURBULENCE NOISE 


09999 FILENAME OUTPUT 

10000 DIMENSION AEV< 46) > AEDC 46) > ZMM< A6> , 2PPC 46> / DCVl C9 1)/DCD1(9 1)> 

10001 & F< 20)> DCV2(9 1) > DCD2(9 1) # STH0SRC9 1) > AEVETCC 2>9 1) / 

10002 & AEDETC<2/9 1)> EL1N(25) 

10010 0UTPUT = ''I TOUT” 

10020 BEGIN FILE OUTPUT 
10030 END FILE OUTPUT 
10100 100 F0RMAT<V) 

10200 READ (’'ITINPUT'S 100) LN/ EM A, EMT/ TPR 
10205 READ ( ”1 TI NPUT"« 1 00) LN# NELMAX/ M/ TI 
10210 READ < "I TINPUT’’^ 100) LN, GAM^ RHO / C S> AI / SI GR 
10220 IBEG«1 
10230 I END- 10 

10240 240 IF <IEND .GT. M) I END-M 

10250 READ C "I T1 NPUT", 1 00) LN^ ( FC I ) , I =I BEG, I END) 

10260 IF (I END .EQ. M) G0 T0 300 

10270 IBEG-IBEG'UO 

1 0280 lEND-I 10 

10290 G0 T0 240 

10300 300 lBEG-1 

10310 I END- 1 0 

10320 320 IF (I END .GT. NELMAX) I END-NELMAX 

10330 READ < "I T1 NPUT", 1 00) LN, ( ELI N( I ) , I =I BEG, I END) 

10340 IF <IEND .EQ. NELMAX) G0 T0 500 

10350 lBEG-IBEG+10 

10360 IEND-IEND+10 

10370 G0 T0 320 

10500 500 WRITE (OUTPUT, 510) 

10510 510 F0RMAT(//32H EMA EMT M TI ) 

10520 WRITE (0UTPUT, 530) EMA, EMT,M,TI 
10530 530 F0PMAT(2F9.3,16, F9.2) 

10535 C-CS 

10536 A=AI 

10540 WRITE (0UTPUT, 550) 

10550 550 F0RMAT(//53H GAM RH0 C A SI GR 

10551 & TPR) 

10560 WRITE C0UTPUT,57O) GAM, RH0, C, A, SI GR, TPR 
10570 570 F0RMAT(F9.3,F9.4,F9.O,F9.2,F8.2,F9.2) 

1 0600 DBL= 1 30. + 4. 3429 45*AL0G< . 1O5*RH0*( C*EMA) ♦♦3*TI 

1 0610 PI =3. 141 5926 

10615 TP1=2.*PI 

10620 G10V2 = (GAM- 1 .)/2. 

10630 G10VG=(GAM-1.)/GAM 
10640 EMR-SQRT(EMA**2-*>EMT**2) 

10650 Tl 1-TPR**G10VG-1. 

10660 T12=l .+ 1 ./(G10V2*EMA**2) 

10670 SR1MM2=SQRT< 1 .-EMA**2) 

10680 SR1MP2=SQRT< 1 .-EMR**2) 

1069 0 CR-EMA/EMR 
1069 5 CR2=0R*CR 
10700 SR-EMT/EMR 
10705 SR2=SR*SR 
10710 SRCR=SR*CR 
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10715 TR=EMT/EMA 
10720 EMRC«EMP./( 5.*CR) 

10730 CL=CR2*T12*T1 l/< SIGR*SR> 

10740 AA=1.-EMA**2 

10750 A=CL*S1 GR*EMA*EMT/(4.*AA) 

10760 B=»CL*SI GR*SR1MP2/(4.*AA) 

10770 C=CL*SIGR/4. 

10780 ALC=TP1 ♦SR1MR2 
10790 BETC-TPI 
10800 CHIC=BETC*EMT 
10810 DELC=TP1*EMA*EMT 
10840 CDP=PI*SI GP,/ 2. 

10845 CDP2=CDP*CDP 
10850 CAE=l ./<EMA*SR1MM2) 

1 09 00 NINC0 2= 10 
109 10 NINC=2*NINC02 
10920 FNINCaNINC 
10930 DELTH = PI/FNINC 
10940 THETA«“PI/2. 

109 50 1MAX=N1NC-M 
1 1000 D0 1230 I = 1>IMAX 
11010 CTH=C0S< THETA) 

11020 EMAMC»EMA-CTH 
1 1030 EMAPC=EMA+CTH 

1 1040 IF (I .GT. (NINC02+1)) G0 T0 109 0 
1 1050 ZMMO )»EMAMC 
11060 ZPP<I)=EMAPC 
11070 AEVa)=CAE/<l.+EMA*CTH>**2 
11080 AED<I )=CAE/( 1 .-EMAfCTTO^^B 
1 1090 1090 STH = SIN< THETA) 

11100 TERM=SP.1MM2*STH*CR 
11110 DCVl < 1 ) *<EMAPC*SR-TEFM)**2 
11120 DCV2(1)=<EMAPC*SR+TERM)**2 
11130 DCDKI )=<EMAMC+SR-TERM)**2 
11140 DCD2CI)=<EMAMC*SR+TERM)**2 
11150 STH0SP< I ) =S7H/SRIMM2 
1 1160 INDEX=I 

11170 IF (I .GT. (NINC02+D) I NDEX = I MAX+ 1 - 1 
11180 AEVETC( l/I)aAEV(INDEX)*CDP2*DCVl(l) 

1 1190 AEVETC<2W)=AEV<INDEX)+CDP2*DCV2CI) 

1 1200 AEDETCC 1 / I ) =AEDU NDEX) *CDP2*DCD1 (I ) 

11210 AEDETC<2/I)=AEDCINDEX)*CDP2*DCD2(I) 

1 1220 THETA=THETA+DELTH 
1 1230 1230 CONTINUE 
11250 00 3010 NEL=UNELMAX 
11255 EL=ELINCNEL) 

11260 WRI TE (0UTPUT. 1261) EL 

11261 1261 F0RMAT<//11H ***** EL= 1 P 1 E 1 2 . 3^ 6H *****) 
1 1270 CPH1«1,./(4.*PI*EL) 

11275 ELC=1 ./<EL*EL) 

1 1300 00 3000 J= UM 

11310 WRITE <0UTPUT, 1320) F<J) 

11320 1320 F0RMATC//7H F(J)*E12.3) 



1 1330 CHI=F<J)*CHI C 
1 1335 CHI 2=CHI*CH1 
I 13i»P SDV*0. 

1 1350 SDD-0. 

1 1360 SQV=0. 

1 1370 SQD=0. 

11380 NVAL=8 
1 1390 NVALP1=NVAL+ 1 
1 1500 D0 2800 NN=»i,NVALP1 
11510 N=NN- 1 

1 1520 IF <N> 1530, 1550> 1 570 

1 1530 1530 SGN»- 1 . 

11540 G0 T0 158 0 
1 1550 1550 SGN=«0. 

1 1560 G0 T0 158 0 
1 1570 1570 SGN=1 . 

1 158 0 1 58 0 EN»N 
11590 AL=ABS(EN>*ALC 
11595 AL2»AL*AL 
1 1600 ALSGNB=AL*SGN+B 
1 1605 BET»EN*BETC 
11610 DEL=DELC*EN 
1 1630 SUMDV=0. 

1 1640 SUMDD=0. 

1 1650 SUMQV=0. 

1 1660 SUMQD=0. 

1 1670 THETAs-PI/2, 

1 1800 D0 2600 I=l,IMAX 
1 1805 AKX=BET*TP.-CHI/EMA 
11810 AKY = -BET-t-CHl*STH0SR(I > 

1 1820 AKY2=AKY^AKY 
1 1830 INDEX = I 

1 1840 IF <I .GT. <NINC02+1>> I NDEX = I MAX+ 1 - 1 

11850 EM=ZMM< INDEX) *081 

11860 ZP=ZPP<INDEX)*CHI 

1 1900 D0 2310 Ll = l,2 

11910 AKX2=AKX*AKX 

1 1920 ELAKX2=ELC+AKX2 

11930 CPHI 1=CPHI/<ELAKX2+AKY2)**2.5 

1 1940 PHIXX=><4.*AKY2+ELAKX2)*CPHI 1 

11950 PHIXY=-3.*AKY*AKX*CPHI 1 

1 1960 PHIYY = < 4.*AKX2+ELC-»-AKY2)*CPHI 1 

1 1965 TPHIXY=2.*PHIXY 

11970 0MR=SI GR*< AKX*CR+AKY*SP)/2. 

1 1980 SRF=1 ./< I. + TPI*ABS<0MR) > , 

1 1990 PHI T=PHIXX%SR2-TPHIXY*SRCR+PHIYY*CR2 

11995 PHI TS*PHI T*SRF 

12000 IF <L1 .EQ. 2) G0 T0 2040 

12010 FNDVP=AEVETC< 1,I)*PHITS 

12020 FNDDP*AEDETC< 1, I )*PHI TS 

12030 G0 T0 2060 

12040 2040 FNDVM=AEVETC( 2, I )*PH1 TS 
12050 FNDDM=AEDETC( 2, I )*PHI TS 
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12060 2060 ZTERM=ZF 

1 2070 AKXA=AKX*AA 

12080 D0 2200 L2=N2 

12090 ZDELAK=ZTERM-DEL-AKXA 

12100 DEN=( AL2+ZDELAK»*2)»*2 

12110 PART= ZTERM* C Z DEL AK* A- AL SGNB) + C*Z DEL AK* CH I ♦ STH0 SR( I ) ♦ AA 

12120 GX=ZTERM*PART 

12130 GY=CHI*STH0SR<I)*PART+AA 

12140 T0=( GX**2*PHIXX+GY+*2*PHIYY+GX*GY*TPHIXY)/DEN 
12150 IF (L2 .EQ. 2) G0 T0 2180 
12160 FNQV= TO* AEV< INDEX) 

12170 G0 T0 2190 

12180 2180 FNQD»TQ*AED< INDEX) 

12190 2190 ZTERM=»ZM 

12200 2200 CONTINUE 

12210 IF (LI .EG. 2) G0 T0 2250 

12220 FNQVP*FNQV 

12230 FNQDP=FNQD 

12232 AKX=BET*TR+CHI/EMA 

12234 ZP=-ZP 

12236 ZM=-ZM 

12240 G0 T0 2310 

1 2250 2250 FNQVM=»FNQV 

1 2260 FNQDM=FNQD 

1 2310 2310 CONTINUE 

12400 FDV=FNDVP+FNDUM 

12410 FDD=FNDDP+FNDDM 

•1 2420 FQV=FNQVP+FNQVM 

12430 FQD=FNODP*FNQDM 

12440 IF ((I .NE. 1) .AND. (I .NE. IMAX)) G0 T0 2490 
12450 FDV=FDV/2. 

12460 FDD=FDD/2. 

1 2470 FQV=FQV/2. 

12480 rQD=FGD/2. 

1249 0 249 0 SUMDV= SUMDV* FDV 

12500 SUMDD=SUMDD+FDD 

12510 SUMQV=SUMQV+FQV 

12520 SUMQD=SUMQD+FQD 

12530 THETA=THETA•^DELTH 

12600 2600 CONTINUE 

12610 SNDV=CHI 2*SUMDV*DELTH 

12620 SNDD=CHI 2*SUMDD*DELTH 

1 2630 SNQV=SUMQV*DELTH 

12640 SNQD=SUMQD*DELTH 

12650 IF (N .GT. 0) G0 T0 2750 

12700 SDV=SDV+SNDV 

12710 SDD=SDD+SNDD 

12720 SQV=SOV+SNQV 

12730 SGD=SGD+SNQD 

12740 G0 T0 28 00 

12750 2750 SDV=»SDV+ 2 . * SNDV 

12760 SDD=SDD•^2.*SNDD 

12770 SGV=SQV+2.*SNQV 



12780 SQD»SQD+2.*SNQD 
12800 2800 CONTINUE 

12802 SDV»SDV«EMRC 

12803 SDD»SDD*EMRC 

12804 SQV»SQV*EMRC 

12805 SQD»SQD*EMRC 

12810 SDVDBR«>4.34294S*AL0G(SDV) 

12820 SDDDBR-4.3429454<AL0G(SDD> 

12830 SQVDBRo4.342945*AL0GCSQV) 

12840 SODDBR«4.3429 4S*AL0G(SQD) 

12850 PV=SDV+SQV 
12860 PD»SDD+SQD 

12870 PVDB=4.342945*AL06<PV)^DBL 
12880 PDDB-4.3429 4 5«AL0G(PD)4-OBL 
129 0 0 WRITE (0UTPUTj29 10) 

129 10 2910 F0RMAT<//4OH SOVDBR SDDDBR SQVDBR 
1 2920 WRITE (OUTPUT^ 2930) SDVDBR# SDDDBR, SQVDBR# SQDDBR 
129 30 29 30 F0RMAT( 4F1 0. 1 ) 

129 40 WRITE C OUTPUT, 29 50) 

129 5 0 29 50 FORMAT! //20H PVDB PDDB) 

1 29 60 WRI TE < 0 UTPUT, 29 70) PVDB, PDDB 
129 7 0 29 7 0 FORMAT! 2F1 0. 1 ) 

13000 3000 CONTINUE 
13010 3010 CONTINUE 
13100 STOP 
131 10 END 


SQDDBR) 
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INLET DISTORTION NOISE 


09990 FILENAME OUTPUT 

10000 DIMENSION SI GN( 2) / ETA< 2) > TPDS< 50> / DELM( 500) ^ PD( 2) ^ PQ( 2) 

10010 0UTPUT = "ID0UT" 

10020 begin file output 

1 0030 END FILE OUTPUT 
10100 100 FORMAT (V) 

10200 BEAD ("I DINPUT'V 100) LN, TPD> I S^ EMA> EMT^ NDP 
10210 BEAD ("IDINPUT"/ 100) LN^ GAM> RHO^ 0/ A/ N> 1 B, SI GR/ TPR 
10220 I BEG-2 
10230 I END- 1 1 

102A0 240 IF (I END .GT. CNDP+D) IEND-NDP+1 

102 50 BEAD ( "1 DI NPUT”/ 1 00) LN, C TPDS< 1 ) / I =I BEG> I END) 

10260 IF (lEND .,EQ. < NDP+ 1 ) ) GO TO 500 

10270 IBEG-IBEG+10 

10280 lEND-IEND+10 

1 0290 GO TO 240 

10500 500 S-IS 

10510 EN-N 

10520 B=IB 

1 0600 WRITE (0UTPUT^61O) 

10610 610 F0RMAT<///36H TPD IS EMA EMT) 

1 0620 WRITE (OUTPUT, 630) TPD> I S, EMA, EMT 
10630 630 F0RMAT(F10.3,I7,2F10.3) 

1 0640 WRITE (OUTPUT, 650) 

10650 650 F0RMAT(//65H GAM PH0 C A N IB 

10651 & SIGR TPR) 

10660 WRITE (OUTPUT, 665) GAM, RHO , C, A, N, I B, SI GR, TPR 
10665 665 F0PMAT(F9.3^F9.4,F9.0, F9.2,2I6,F9.3, F9.2) 

10670 PI*3. 1415926 

10675 IMAX-500 

10676 I MAXI -I MAX- 1 

10677 FI MAX 1=1 MAXI 
1 0680 JMAX=NDP+2 

1 0685 JMAXl-JMAX- 1 
1068 6 FJMAX1=JMAX1 
1068 7 TPDS(1)=0. 

10688 TPDS( UMAX) -0. 

10690 DELM(1)=0. 

10700 DELM( IMAX) =0. 

10710 G10V2 = ( GAM- 1 .)/2. 

10720 G10VG=( GAM-1 .)/GAM 
10730 Tl= 1 « + G 10V2*EMA**2 
10740 T1 1=TPR**G10VG- 1. 

10750 T12= 1 .+ 1 ./( GIO V2*EMA**2) 

10760 DO 800 I=2,IMAX1 
10762 FI=I 

10764 TERMI =( FI - 1 .)/FIMAXl 
10766 J1=TERMI*FJMAX1+ 1. 

10768 J2=J1+1 
1 0770 FJ 1=J 1 
10772 FJ2-J2 

10774 TERMJ 1»(FJ 1- 1 .)/FJMAXl 
10776 TERMJ2=( FJ2- 1 .)/FJMAXl 
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1 0778 TPDI *TPDS( J 1 ) + < TPDSC J2) -TPDSC J !) )*( TEPMI -TEPMJ 1 ) /( TERM J 2- TERMJ 1) 

10780 TPDD=TPD*TPDI 

10785 T2 = < 1 .-TPDD)**G10VG 

1 079 0 DELMC I ) »EMA-SQPT< < T1*T2- 1 .)/G10U2> 

10800 800 C0NT1NUE 

10810 EMR=SQPT(EMA**2+EMT**2) 

10820 SR1MM2»SQRTC 1 . -EMA*«2) 

10830 SR1MR2=S0RT< 1 .-EMP»*2) 

108A0 IM1N«( 1 .-EMT/SRlMM2)oEN*B/S*i999 
1 0850 IMAX»( 1 . + EMT/SR1MM2)*EN*B/S 
10860 ALSRR*ATAN<EMT/EMA) 

10870 VR=C*EMA/C0S< ALSRR) 

1 088 0 C0NST»< VR*VR*SI Gfv 2. ) **2* A*RH0* 7 A6. / C C* 32 • 2* 550 . ) 

1089 0 SR«S1N< ALSRR) 

109 00 CXI = SR 

10910 TXI =-C0S< ALSRR)/SR 

10920 CLT=C0S( ALSRR)**2*T12*T1 l/< SIGR*SR) 

10930 AA»1 .-EMA**2 
10940 AP*CLT*EMA*EMT/AA 
109 50 BP=CLT*SR1MR2/AA 
109 55 CP-CLT 
10960 BB»EN*B*EMT 
10970 BB2=BB*BB 
1 0980 CC=EMA*BB 
1099 0 DD=EMA/BB 
1 1000 ALPH=EN*B*SR1MR2/AA 
11010 SPUD=0. 

1 1020 SPDD=0. 

1 1030 SPUQ*0. 

1 1040 SPDQ=0. 

1 1050 SIGN< 1)=1 . 

1 1060 SI GN( 2) =- 1 . 

1 1070 C0M»PI*SR*SI GR*S/B 
1 1080 CCLC = PI*SIN(2.fALSRR)/EMA 
11090 WRITE (0UTPUT, 1800) 

11100 D0 1900 I=IMIN*1MAX 
11110 rial 
11120 M»N*I B-I*I S 
11130 EM=M 

1 1200 DEL TH=2.»PI/<S*FI MAXI) 

11230 THETA=DELTH 
1 1240 SUM 1=0. 

1 1250 SUM2=0. 

1 1300 D0 400 K=2/1MAX1 
11310 TERM*FI*S*THETA 
11320 F1=DELM(K)*C0S<TERM) 

11330 F2=DELM(K)*SIN<TEPM) 

1 1340 SUM1 = SUM 1+Fl 
1 1350 SUM2=SUM2+F2 
1 1360 THETA=THETA+DELTH 
1 1400 '400 CONTINL'E 
11410 TERM1=DELTH*S/C 2.*PI ) 

1 1420 Cl R=TEPM 1*SUM 1 
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1 1430 Cl I =TERM 1*SUM2 
11440 Cl =SQRT( Cl R**2+CI I +*2) 

1 1450 DEL = SQRT( BB2-AA*EM**2) 

1 1460 T0PD=(CXI/DEL)*CCLC*C1 
11470 EMTXI=EM*TXI 

1 1480 T0PQ=CI/<< ALPH**2+(DEL/AA)**2)*EMR> 

1 1490 0M=F1 *C0M 
1 1500 APPa. 18 11. 

1 1510 SEARS=SQPT<<APP+0M)/(APP-KP1*APP+1.)*0M+2.*PI*0M*0M)) 

1 1600 D0 1700 K» N 2 

11610 CAY=( SI GN<K)*CC+DEL)/AA 

11620 A0VFDaT0PD*CAY*< l. + EMTXl*SIGN(K>/CAY)*SEARS 
1 1630 Fl=CAY/S0RT(CAY**2-fEM*>f2) 

11640 F2= 1 .+SI GN<K)*CAY*DD 

11650 ETACK>a. 5*<F1-SIGN(K)*EMA)/F2 

1 1660 AGVFQaCAY+C AP«'DEL*CAY/AA-ALPH*BP*CAY*S1 GNCK)+CP*SI GN(K)*DEL* 

11661 & EM/AA)*T0PQ/DEL 

1 1670 PD<K)=ETA(K)*A0VFD**2 

11680 PGKK) =ETA<K) *A0 VFQ**2 

1 1700 1700 CONTINUE 

11710 PUD=PD< 1)*C0NST 

11720 PDD=PD( 2)*C0NST 

1 1730 PUQ = PQ< D+C0NST 

11740 PDQ=PQC 2)*C0NST 

1 1750 PUDDBal3O.+ 4.342945*AL0G(PUD> 

1 1760 PDDDBa 1 3O.+4.342945*AL0GCPDD) 

1 1770 PUQDB=13O.+4.342945*AL0G(PUQ> 

1 1700 PDQDBa 1 3O. + 4.342945t'AL0G<PDQ> 

11800 1800 F0RMAT <//50H I M PUDDB PDDDB PUODB 

11801 A PDQDB) 

118 10 WRITE C0UTPOT/ 13 20) I , M/ PUDDB# PDDDB/ PUQDB/ PDQDB 
11820 1820 F0RMAT (2I5#4F10.1) 

1 1830 SPUD=SPUD+PUD • 

1 1840 SPDDaSPDEH’PDD 
11850 SPUQaSPUQ+PUQ 
I 1860 SPDQ=SPDQ+PDQ 
1 1900 1900 CONTINUE 

11910 SPUDDBal30.+4. 3429 45* AL0G( SPUD) 

11920 SPDDDB= 1 30.*4.342945*AL0G( SPDD) 

1 1930 SPUQDB- 1 30. + 4. 34 29 4 5*AL0G( SPUQ) 

11940 SPDCDB* 1 30.+4. 342945*AL0G( SPDO) 

1 19 50 WRITE (OUTPUT, 19 60) 

11960 1960 FORMAT </41H SPUDDB SPDDDB SPUQDB SPDQDB) 

11970 WRITE (0UTPUT# 1980) SPUDDB# SPDDDB# SPUQDB# SPDQDB 

1 1980 1980 FORMAT (4F10.1) 

2 0000 STOP 
2 0010 END 



APPENDIX 2 


LIST OF SYMBOLS 



B 

C 


d 

j 

L/D . 

M M^ and M. 
a r t 

P 

P02'^Pq1 

Pt 

R 

S 

^or ^02’ ''^1 



speed of sound. 

sometimes used to denote radius of fan annulus being 
studied. 

number of blades in the rotor, 
lift coefficient of the rotor, 
rotor chord. 

specific heat at constant pressure, 
coefficients of inlet distortion, 
pitch of rotor. 

integral length scale of turbulence/rotor pitch, 
axial, relative and wheel tip Mach numbers, 
pressure. 

pressure ratio across rotor 
local total pressure 
gas constant 

number of lobes in distortion. 

total temperature upstream and downstream of rotor, 
temperature ahead of rotor. 

axial velocity entering rotor. 

wheel velocity. 

relative velocity through rotor, 
stagger angle of rotor. 

Kronecker delta, 
delta function, 
density, 
rotor solidity. 
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FIGURE 1 


CALCULATION OF ROTOR POTENTIAL FLOW FIELD 
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FIGURE 3. INLET TURBULENCE - ROTOR NOISE 










0.55, M. = 0.825, Rest as- in Figure 3(a) 



FIGURE 4 . INLET TURBULENCE - ROTOR NOISE (Contd. 
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0.35, M. = 0.7, Rest as in 4(a) 







FIGURE 5. EFFECT ,0F TIP SPEED AT CONSTANT PRESSURE RATIO 



0.5, Rest as in 5(a) 












FIGURE 6. INLET DISTORTION NOISE ■ FIGURE 7. INLET DISTORTION NOISE 



QUADRUPOLE 





o 80% Speed 
A 70% Speed 


Data from NASA 
TMX-2191 


y 60% Speed 
O 50% Speed 
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DIPOLE NOISE DUE TO INLET DISTORTION 
(THEORETICAL) 

QUADRUPOLE NOISE DUE TO INLET DISTORTION 
(THEORETICAL) 







FIGURE 9 . COMPARISON WITH (5); ROTOR 2 




one lobe 



V^^Pt^max = TPDS 


FIGURE A1 

DEFINITION SKETCH FOR INPUT FUNCTION TPDS 
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